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In many numerical implementations of the Cauchy formulation of Einstein's field equations one 
encounters artificial boundaries which raises the issue of specifying boundary conditions. Such 
conditions have to be chosen carefully. In particular, they should be compatible with the constraints, 
' yield a well posed initial-boundary value formulation and incorporate some physically desirable 

, properties like, for instance, minimizing reflections of gravitational radiation. 

■ Motivated by the problem in General Relativity, we analyze a model problem, consisting of 

04 ' a formulation of Maxwell's equations on a spatially compact region of spacetime with timelike 

boundaries. The form in which the equations are written is such that their structure is very similar 
^ . to the Einstein-Christoffel symmetric hyperbolic formulations of Einstein's field equations. For this 

' model problem, we specify a family of Sommerfeld-type constraint-preserving boundary conditions 

and show that the resulting initial-boundary value formulations are well posed. We expect that 
these results can be generalized to the Einstein-Christoffel formulations of General Relativity, at 
least in the case of linearizations about a stationary background. 
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For a number of years there has been an effort in the field of General Relativity to obtain a set of boundary conditions 
which are consistent with constraint propagation, physically reasonable and stable in the sense that they yield a well 
posed initial-boundary value formulation. Such a set is essential when integrating Einstein's field equations on a 
domain with artificial timelike boundaries on a time scale of astrophysical relevance, as is the case in most numerical 
^ i' simulations of the binary black holeproblem. Besides the success of 0] for the case of their frame formulation, and 
1^' in spite of a lot of efforts Q| no boundary conditions which satisfy all the above properties 

^ • have been found for the more commonly used tensorial formulations. In an attempt to understand this problem we 
consider in section ^1 a toy model which in many aspects resembles some of the evolution systems currently employed 
for numerical evolution in General Relativity. This model is just Maxwell's equations in a 3 -I- 1 decomposition, but 
written as a first order system in terms of the vector potential, the electric field, and all first order spatial derivatives 
of the vector potential. The allowance of all first order spatial derivatives of the vector potential as variables, and 
' not just the antisymmetric, gauge-independent ones which describe the magnetic field, makes the situation similar 
to the one on the above mentioned systems, where all derivatives of the metric are indiscriminately promoted to 
evolution variables, regardless of their gauge dependence. Thus, the Maxwell system, which in gauge-independent 
variables (electric, and magnetic fields) has a nice symmetric hyperbolic initial value formulation as a set of evolution 
equations, here acquires most of the pathologies one encounters for the general relativistic systems. In particular, 
unphysical constraints appear, and the constraints propagate with nontrivial speeds at timelike boundaries. The 
indeterminacy of the evolution system under the addition of linear combination of constraints is used to construct a 
two-parameter family of evolution systems. For some values of the parameters the system is symmetric hyperbolic, 
for some others it is only strongly hyperbolic, while still for others it is neither, and so ill posed. This system has 
been used by a number of researchers in the past, in particular see 

The construction of boundary conditions for this toy model is easier than for the case of Einstein's equations, and 
we can look at it in full generality. In particular, in the toy model, we know without ambiguities what the physically 
relevant boundary conditions are which lead to non-incoming radiation conditions, and so we can impose them, 
together with some conditions ensuring constraint propagation. As shown in section UTTI for a range of parameters this 
set of boundary conditions satisfies the so called Kreiss condition |il2t. il"3L IT^ and so they yield a well posed problem 
for trivial initial data. But it turns out that, surprisingly, this set of conditions does not yield well posed problems in 
the traditional sense, that is when non-trivial initial data is allowed. This is shown by a counterexample consisting of 
a solution to the equations which docs not satisfy the expected energy estimates. This counterexample also indicates 
the source of the problem which is, basically, the existence of static solutions which are represented by linearly in time 
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growing gauge potentials, and motivates a new gauge condition. By using this new condition, we are able to show 
well posedness in a Hilbert space that controls the norm of the main variables and the constraint variables. In 
section Hvl we derive a priori estimates. The main idea is to start with an estimate for the gauge-invariant quantities 
instead of estimating the norm given by the symmetrizer of the evolution system, which is gauge dependent. Based on 
these estimates we show existence of solutions in section^ Since the new gauge condition is imposed by some elliptic 
equations, some elliptic theory is needed. For the proof of existence we employ the abstract theory of semigroups, 
for which the elliptic aspects of the problem are taken care of in a natural way. In section IVII we summarize our 
results and discuss some of their implications for the construction of boundary conditions for tensorial formulations 
of Einstein's field equations. We have also provided an appendix with the details of the elliptic theory needed; the 
material there is standard. The main results and the counterexample that motivates the new gauge condition are 
stated in the next section. 



II. MODEL PROBLEM AND MAIN RESULTS 



In this section we present our model problem which consists of a system of evolution equations with constraints. We 
start with the case without boundaries and state in which sense the resulting initial-value problem is well posed. In 
particular, we introduce a Hilbert space Ti which controls the norm of the main and constraint variables. We then 
consider the presence of artificial boundaries, discuss boundary conditions and show by an explicit counterexample 
that the resulting initial-boundary value problem is not well posed in the expected space Ti. This motivates a new 
gauge choice. Finally, we state the main result of this article which proves well posedness in Ti, of the initial-boundary 
value problem for this new gauge choice. 

Let ri = M'^, and denote by V the covariant derivative with respect to the Eulerian metric h = Sijdx^dx^ on R^[3^. 
We are interested in the following evolution system on fl: 

dtp = VV, , (1) 
dtA^ = E,+V^(l), (2) 
dtEj = V'Wij - (1 + a)\7^Wji + aVjW + Jj , (3) 



where p and are scalars on Xl, ^4^, Ei and Ji are one-forms on f2, and Wij a two-tensor on Q, with trace W = h^^Wij. 
a and /? are two parameters which determine the dynamics off the constraint hypersurface which is defined by 

C = p- V'^Ek = 0, (5) 
C^P ^ W,, - V,A, = 0. (6) 

Physically, the system 12l3l4l5lt)|l is equivalent to Maxwell's equations on Minkowski space, where represents the 
electrostatic potential, Ai the vector potential, Ei the electric field, (Bi) = {W23 — W32, W31 — M^ia, W12 — W21) the 
magnetic field, and p and Jj represent the charge and current density, respectively, which obey the continuity equation 
(Q. Notice that for given current density, Eq. can be integrated separately in order to obtain p, which can then 
be used as a source function in order to integrate the equations (|2I3I4|I . However, we will find it more convenient 
to interpret p as a field being evolved along with the fields Ai, Ei, Wij since in this case the constraint variable C 
depends linearly on the evolution fields. We assume that the electrostatic potential and the current density are a 
priori given. The motivation for introducing the fields Wij, which represent all the first order spatial derivatives of 
the vector potential, instead of using the magnetic field is to obtain a system of equations whose structure is similar 
to the one of the Einstein-Christoffel formulation of Einstein's field equations [iSl iWL IitI Il8| . In the latter, one has 
evolution equations for the components of the three-metric gij , the extrinsic curvature Kij , and some symbols dkij 
that are linear combinations of the Christoffel symbols. These fields are subject to the Hamiltonian and momentum 
constraints, and to the constraints dkij — dkQij = 0. The structure of the equations is very similar to the ones in our 
model problem with the correspondence gij Ai, Kij ^ Ei, dkij ^ ^iji ^ 4'i where TV* is the shift vector field. 

In the following, we impose the condition a ■ (3 > which is a necessary and sufficient condition for the evolution 
system (I1I2I3I4II to be strongly hyperbolic which in turn yields well posedness of the associated Cauchy problem in 
L^(ri) 13] . The constraints' propagation is described by the following evolution system 

dtC - -aV'^Ck , (8) 
dtCk = -l3VkC, (9) 
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which is a consequence of the evolution system (|1I2I3I4() . Here, we have introduced the new constraint variable 
Cfc = 2/i*W[fcC|j^-' = VfeVF — V^Wkj in order to obtain a first order evolution system. A simple energy estimate 
shows that the unique solution to Eqs. (|7I8I9|I with trivial initial data is trivial, and therefore, the constraints remain 
satisfied if satisfied initially. Summarizing, we have the following 

Theorem 1 (Well posedness of the Cauchy problem) For a/S > the constrained Cauchy problem associated 
with < 1^516]) on rt = MP is well posed in the following sense: Let t > be an arbitrary fixed constant with di- 

mension of length, and write u(t) — {Tp{t),T~^Aj{t), Ei{t),Wij(t)) (main variables), v{t) = {T~^C^^\t),C(t),Ck{t)) 
(constraint variables) and j(t) = {Ji{t),T'SI^ Ji{t),T^^'S/ i4>,'S/ iS/ jcj)) (source functions). Given smooth source functions 
j{t) G i^(ri), t > such that 1 1~> j{t) is continuous, and given smooth initial data uq G L'^{n), there exists a unique 
solution u{t) G L'^{n) of the evolution system mijl.'jl il) with u(Q) = uq. Furthermore, the solution obeys the estimate 



"o||l2(o) +t ||j('S)llL2(n)(is 



(10) 



for some constants a, b. 

If the initial data is such that vq G L^(il), then v{t) G L'^{fl) and we also have the estimate 

Mt)\\hin)<ce''/-\\vo\\h^a), (11) 

for some constants c, d. In particular, this implies that initial data which satisfies the constraints initially automatically 
satisfies the constraints at later time. 

Remarks: 

1. The constant r is an artificial length scale that is introduced in order to avoid adding quantities which have 
different units. By choosing this scale to be arbitrarily large we can make the growth rate 1/r in the estimates 
as small as we like. 

2. The precise sense in which there exists a solution u{t) G L^(ri) is in the sense of a strongly continuous semigroup 
in L^(ri) whose generator is determined by the operator on the right-hand side of 



3. In the following, we will replace the Hilbert space {u — {rp, r ^Ai,Ei, Wij) G L^{n)} by the Hilbert space 

H = {u= {tp,t-^A,,E,,W,,) G L^{n) : v = {r-^ciP ,C,Ck) G L\n)} (12) 

with scalar product (ui , U2)h — (ui, U2)i,2(si) + r'^ivi, V2)L^{n)- This Hilbert space controls the L^ norm of the 
main variables and the constraint variables. In this case, we replace the estimate ()l(Jf) by the estimate 



Wum'^Kae"'/^ 



\\uq\\h + t / \\j{s)\\l2(^a)ds 







(13) 



The space Ti. will be important when boundaries are present, where we will be able to show well posedness in 
n (see Theorem 12 below) but not in L^. That is, we will be able to derive an estimate which involves the L 
norms of the main and the constraint variables, but not the main variables alone. 

The purpose of this article is to analyze the constrained evolution system (|1I2I3I4II . (|5I()|I on a bounded domain 
of R'^ with C°° boundary dil. In the following, n — nidx^ denotes the outward unit one-form to dVl. We also 
introduce the projection operator = — UiU^ on the tangent space of dfl. For technical reasons, we assume that 
the extrinsic curvature of the boundary surface, defined by Kij — H[ H^V rris^ is positive semi-definite at each point 
of the boundary. 

We wish to specify boundary conditions on dfl which imply a well posed initial-boundary value problem in Ti.. 
In particular, these conditions have to be specified such that the constraints are propagated which means that the 
constraints (|5I6|I should hold at later times if satisfied initially. We shall call boundary conditions which have this 
property constraint preserving. In order to specify such conditions we impose homogeneous maximally dissipative 
boundary conditions for the constraints' propagation system (|7I8I9|) . which means that we couple the in- to the 
outgoing characteristic fields with a smooth coupling function c on dfl with the property that |c| < 1: 



(14) 
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where 

P P 

are the in- and outgoing fields, respectively. A simple energy argument (see the next section) shows that the unique 
solution to Eqs. H7I8I9|) with the boundary condition l|14l) and trivial initial data is trivial. Therefore, the boundary 
conditions H14() are constraint preserving. Furthermore, we impose the boundary conditions 

= 5, 'w^^ + 5, , w'^^^ = 'Ej ± n''{Wk^ ~ W,,.), (15) 

where S'^ is a smooth matrix-valued function on dVl with the properties that ■'n' = 0, S^ ■'rij — and S.^ ^ ^ < 
H^'^ if gi = and H^^ ^^j'^ — ^H"^^ for some 5 < 1 otherwise. The function gi is a smooth function on with 
the property that rfgi — 0. Physically, the boundary condition H15() permits to control the normal component of the 
Poynting vector P„ = —n'^E^{Wij — Wji): When gi — 0, this condition makes sure that P„ is nonnegative, meaning 
that the total energy flux through the boundary is nonnegative. When S'j =0 the data gj permits to introduce a 
wave which travels in normal direction towards the boundary. 

With these boundary conditions one would expect to have an estimate, and the corresponding theorem of uniqueness 
and existence, of the following form: 




for some constants a, b. But this expectation is false, as shows the following counterexample. 

A. Ill posedness in 

Suppose that the source functions (f> and Jj vanish identically, and consider the following family of solutions: 

p = 0, 

= V,/, (16) 

where / is a smooth, time- independent, harmonic function (V'^Vfe/ — 0). This family satisfies the constraints (|5lt)|l . 
the evolution equations H1I2I3I4() and the boundary condition H14() . It has the initial data 

and satisfies the boundary conditions H15() with boundary data given by 

This family of solutions does not obey the desired estimate in Ti since Wtj depends on second derivatives of / whereas 
the initial and boundary data depend only on first derivatives of /js^- Physically, the solution (|16|l represents an 
electrostatic solution in a "bad" gauge corresponding to a boundary with charge density proportional to n'Vi/. This 
gauge is "bad" in the sense that it is not adapted to electrostatic solutions since it requires the electric potential to 
be zero. As a consequence, the potential grows linearly in time. Notice that the solution (|16l) is trivial in the absence 
of boundaries since (if one requires the initial data to decay sufficiently fast at infinity) / has to be zero. This is 
compatible with the fact that in those cases one can show well posedness in since the evolution system is strongly 
hyperbolic (see Theorem^. 

B. Main result 

The above counterexample shows that the simple gauge condition (/) = does not lead to a problem that is well 
posed in TL. This counterexample motivates the following condition on 0: 



V'^Vfc^ = -V^Ek , n^Vk<t> = -n^Ek on dQ.. 



(17) 
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This condition precludes the type of solution leading to the counterexample, and, as we will state shortly, allows for 
the formulation of a well posed initial-boundary value problem. Notice that the gauge condition IjlTI) guarantees that 
dt{n'^Ak) — on dft. In particular, this implies that the solution satisfies n'^Ak = on dfl if the initial data satisfies 
this condition. In order to state the main result of this article, we introduce the Hilbert space 

n' = {uen: n^Ak on dVL}. (18) 

Theorem 2 (Well posedness of the initial-boundary value problem) Consider the constrained evolution sys- 
tem UlSi<^ Jj). 1,516]) on a bounded domain il with C°° boundary and boundary conditions {l^ll,^) . where if) is determined 
by the elliptic system j j7| ). Assume that a(3 > Q, and that the extrinsic curvature of the boundary surface, Kij, is 
positive semi-definite at each point of the boundary. Assume that the matrix-valued function S- appearing in the 
boundary condition J^15}) has the form S^'' = sH^-' for some function s G C°^{dfl), and assume that \c\ < 1, \s\ < 1 
on dn. 

Then, given smooth source functions j{t) = {Ji{t),W^Ji{t)), t > 0, in L'^{il), such that t t-^ j{t) is continuous, and 
given smooth initial data uq G Ti' and smooth boundary data gi{t) in L^{dfl) with appropriate compatibility conditions 
(see Eqs. Hyil9<l}} below) at {t = 0} x dfl, there exists a unique solution u(t) £ Ti.' of i) ^14-115]) with 
u(0) = uq. Furthermore, the solution obeys the estimates 



Mh + ^ I / ||g(s)iii 







(19) 



vm'Un)<<^e'"/^\\vo\\Un)- (20) 



Remarks: 



1. Again, the precise sense in which there exists a solution u{t) £ Ti.' is in the sense of the existence of a strongly 
continuous semigroup in Ti,' . This is discussed in section Ivl 

2. The assumptions on the nonnegativeness of and on the special form of the coupling matrix can be 
weakened, but the proof is technically more complicated in those cases. 

In the next section, we start with some preliminary investigations of the initial-boundary problem (|1I2I3I4|) , (|14I15I) 
which are valid in the "high-frequency limit" . That is, we consider the case where 57 is a half plane and derive some 
necessary conditions for well posedness, using Laplace and Fourier transformation techniques. A complete proof of 
Theorem 13 is given in the subsequent sections. 

III. CONSTRAINT-PRESERVING BOUNDARY CONDITIONS 

In this section we construct a family of constraint-preserving boundary conditions for the case where = {(x, y, z) G 
R'^ : a; > 0} is the half space. We apply Laplace-Fourier techniques in order to derive necessary conditions for 
well posedness of the resulting initial-boundary value problem. We also derive a stronger condition, known as the 
Kreiss condition, which yields an estimate in the case of trivial initial data and trivial sources. However, the 
counterexample presented in the previous section shows that this result cannot be generalized to the case of non- 
trivial initial data when the temporal gauge = is adopted. This shows that the verification of the Kreiss condition 
for constraint-preserving boundary conditions does not necessarily yield well posedness in the space . For simplicity 
we only consider the source-free case here, where p = 0, Jj = 0. We also neglect the evolution equations for Ai and 

Cj'J^'' in this section since they can be integrated separately once a solution of (|3I4|I and (|8I9|) has been obtained. 

In order to construct constraint-preserving boundary conditions, we first find the characteristic speeds and fields 
with respect to the unit outward normal (ui) — (—1,0,0). For the evolution system H3I4|I these are 

±r, vi^^ ^ T ^ W\ (21) 

±1, v'-^^ =EaT {W^a - (1 + a)WA.)) (22) 
0, Wax (23) 

0, xab = Wab~\5abW% (24) 

0, K = l3Wxx-{l + l3/2)w\, (25) 
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where here and in the following, Capital indices refer to the directions y and z which are transversal to the boundary, 
5ab = \ '\i a — B and zero otherwise, and r = yja(3 > 0. For the constraints' propagation system (jHEJ, the 
characteristic speeds and fields are 



±r, -l/(±) = C±^a (26) 



r 

0, Ca- (27) 



Defining the energy norm 



taking a time derivative, using Eqs. (|8I9|1 and integration by parts, we find 



= J / fc^' + ^ C'^Ck ) d^x, (28) 



= -a f {CV'^Ck + C'VkC) d'x = a ( CC, dydz {(V'^+^f - {V^-'^f) dydz. 



(29) 



Therefore, if we impose the boundary condition 

= c y(-) (30) 

with |c| < 1, it follows that the positive definite quantity fc cannot increase; in particular a solution with initial data 
such that C = 0, C/c = satisfies C = and Cfc = at later times as well. Re-expressing condition (|3U|I in terms of 
the variables of the main system yields 



= r(F(+^ -cV^"^) 
= (9t(4+^ +0 1;^")) + V 



U\ - c)Ea + (1 + c)F^A - ^(1 - c)Va<; 



(31) 



where -Fl^yi = W^a — Wax, and where we have used the main evolution equations (|3I4() in order to trade normal 
derivatives (9j.) by derivatives that are tangential to the boundary {dt, Oa)- We can interpret this equation as an 
evolution equation for vi'^^ + c vi at the boundary. Since the main evolution system has three ingoing modes, we 
need two more boundary conditions. One possibility is to require a maximally dissipative boundary condition on the 
transversal fields 

v^+^ ^ Sj" v^B^ + gA , (32) 

where gy, gz are a priori specified functions on the boundary and where the 2x2 matrix S — (S^^) satisfies S'' S < 1. 
From a physical point of view it is more convenient to require 

Ea - F^A = SI{Eb + F,B)+gA, (33) 

since this controls the radiation flux through the boundary, as discussed in the previous section. In order to explore 
both possibilities, we shall analyze the family of boundary conditions consisting of Eq. (|30() and 

[Ea - WxA + (1 + a)WAx] - [Eg + W^b - (1 + a)WBx\ + gA , (34) 
where the parameter a is equal to either a or zero. 

A. Particular cases 

Here, we assume that the parameters a and (3 satisfy — 2<q;<0, /?<— 2/3 which makes sure that the evolution 
system H3I4|) is symmetric hyperbolic. Suppose that </) satisfies V^V a4> = at the boundary. If we choose c = — 1 
and — ~^A^ boundary conditions reduce to 

4+^ = 4"^ + 5, v'^A^ = -v^A ^ + gA , (35) 

where, given gA, the source function g is obtained by integrating the equation dtg = —rP~^'V^gA at the boundary. 
Once g is determined the resulting initial-boundary value problem is well posed since we specify maximally dissipative 
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boundary conditions for a symmetric hyperbolic system. Notice that v)^ + v)[ = 2Ea, so if gA = this means that 
the boundary is a conductor. 

If we choose, instead, c = 1, = and a = a, we obtain 

vi+^ - -vi-^ + 9, v^^^ = v^A^ + gA , (36) 
where now g is determined by the following evolution system at the boundary (given gA and V^^ at the boundary): 

dtg + 2aV^WA. - V V , (37) 
StW^Ax- ^Va9 - VaV,0. (38) 

Since — 2 < a < 0, this boundary evolution system is well posed, and it can be integrated separately to obtain the 
boundary function g (and the zero speed field Wax)- The functions g and gA are then used as boundary data to 
integrate the bulk system, and the resulting initial-boundary value problem is well posed. However, in this case, the 
physical interpretation is less clear since the free data gA = —2FxA + ^uWax is not gauge- invariant. 

These two sets of well posed constraint-preserving boundary conditions have been generalized to the linearized 
Einstein-Christoffel formulation of Einstein's equations 0. However, in order to construct a radiative- type boundary 
condition, as described in the previous section, we need to choose Sa = and a = in Eq. H34|) . and in this case the 
well posedness of the resulting system is more involved since the resulting boundary conditions are not in maximal 
dissipative form. We analyze this in the next subsection. 



B. Fourier-Laplace analysis 



Since the equations and boundary conditions are linear and have constant coefhcients, we can solve the initial- 
boundary value problem by performing a Laplace transformation in time and a Fourier transformation in the spatial 
directions which are tangential to the boundary. In other words, we are considering solutions of the form /(a;)e''*"'"'"-*^ , 
where s is a complex number with positive real part, {to a) = (!^v,^z) is a real two- vector and / € L^(M+). For given 
boundary data there must be a unique such solution for each Re(s) > and ojai otherwise the system admits modes 
that grow like e"* where Re(s) can be arbitrarily large, and the system is ill posed [Tslfl^. 

Performing the Laplace-Fourier transformation, and assuming trivial initial data and = 0, we obtain the following 
system of ordinary differential equations 



sEb 

sWxx 
sWxA 

sWax 
sWab 



adxW\ -f iw^ [Wax - (1 + a)WxA] , 

dxQB + iLo'^WAB - «(1 + a)Lo'^WBA + iaujB{Wxx + W\), 
{l + f]/2)dxEx + '-^uj^EA , 

dxEA , 
iujaEx , 

iLUAEB + ^ Sab {dxEx + ito^ Ec) , 



(39) 
(40) 

(41) 

(42) 
(43) 

(44) 



where we have introduced Qa — WxA ~ (1 ^ 
their evolution equation becomes algebraic: 



ct)WAx- From this we can eliminate the variables with zero speed since 



Wax 

XAB 



i^x , 

S 



loaEb - - Sab^^^Ec 



2s 



(2 + 3(3)uj^Ea . 



(45) 
(46) 
(47) 



Suppose \uj\ — ^ S'^^LOA^B 7^ 0. Let Qja — '^aI\^ 

Introduce the rescaled variables C = s/|a;| and ^ = l^^l^;, and define E^ = Ea(^^, E^j — Eat) 
Qri = QaV^- In terms of these variables, the evolution system decouples into the following two blocks 



and let fiA be a unit two- vector which is orthogonal to loa- 
and define E,., = Ea(^^. Er, = EaV^. Quj = Qai^"^, 



C\( E,, 



(48) 
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c + 







^(l 



C 

- a) 



V 













z(l 



/=) 



g 2/3+1 




(49) 



where Q = aW^j^. The first matrix has the eigenvalues ±A, where A — y/C'^ + 1. The second matrix has the 
eigenvalues ±A and ±fi, where ^ = ^/(^/{aJ3) + l (the sign of the square root is chosen such that for Re(C) > 0, 
Re(A) > and Rc(/i) > 0). The solutions which are bounded as a; ^ oo are given by 



E,j 

Qri 



-c 

A 



(50) 




c 







(72 



_z_ 



(51) 



where cto, cti, (72 are integration constants. A necessary condition for the system to be well posed is that the boundary 
conditions must determine these constants uniquely. Using the relations (|45I46I47|I and the previous notation the 
boundary conditions H3UI34|I yield 



(1 - c)[iaE^ ~ CQ] + (1 + c)^^ 



c- 



= 0, 



(y. — 0, 

(1 - d)E^ -{l + d)[Q^+ I—— E.,\^g^ 



{I - d)E^ - {I + d)Qr, = , 
; assumed tha 

(|50I51|I we obtain the following equations 



(52) 

(53) 
(54) 



where for simplicity we have assumed that 5*^ = dS^ is diagonal. Plugging into this the general decaying solution 



(55) 



/ d_CA - d+(C^ - a) d^C + ad+fi \ f cri \ _ f -iC9uj \ frR\ 

[ c.c-V^c+f,)[<j,/f,) - [ j' ^^^^ 

where c± — c ± I, d± = d ± 1. In order to analyze in what cases these equations uniquely determine the constants 
ctq, CTi and tT2 we use the following 

Lemma 1 Let P > 0, A, B £ M., {A, B) ^ (0, 0), and consider the function 

^ : {Re{0 > 0} ^ C, V^(C) = A C - B^C^+P\ 

where we choose the branch such that Re{^J + P^) > for Re{C) > 0. 

Then, ■0 has zeroes if and only ifA>B>OorA<B<0. Furthermore, \ip\ is uniformly bounded away from zero 
if and only if A ■ B < 0. 

Proof: If S = the Lemma is trivial. So let _B 7^ and rescale C, 4' and A such that V'(C) — ^ ^^C^ + 1. Suppose 
= for some Re(C) > 0. Then it follows that {A^ - 1)C = 1 and thus > 1. Therefore, C = l/V^^ - 1 and 
a/C^TT = \A\/VA'2 - 1. This satisfies V(C) = if and only if A > 1 is positive. This proves the first statement of 
the Lemma. 

In order to prove the second assertion of the Lemma we restrict ourselves to the case A < 1 since otherwise i/j has 
zeroes. Notice first that for large \(\, tpiC) = + hence is not bounded away from zero if A = 1. So 

let A < 1, and let (j be a sequence with Re(Cj) > such that ■0(0) ~^ 0- This sequence is bounded for j large enough 
since otherwise A = 1. Therefore, there is a subsequence that converges to some C* S C. We must have Re(C*) — 
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since otherwise C* is a zero of -0. So we must determine in what cases zero lies on the boundary of the image of V-" in 
C. For C = i cos a, < a < tt, we have vCM-l = sina. Next, let 7 > 0, e > and C — i cosh 7 + e. Then 



+ 1 = y- sinh^ 7 + 2i cosh(7)e + = +i sinh7 + coth(7)£ + O(e^), 



where we have chosen the sign such that Ke{y^(^ + 1) > for e small enough. Therefore, the boundary of the image 
of ip can be parametrized by 

ip(icosa) = — sina + lAcosa, < a < tt, (57) 
'0(±icosh7) = ±i(Acosh7 — sinh7), 7 > 0. (58) 

We see that the boundary contains zero only if < ^ < 1- Therefore, \ip\ is bounded away from zero if and only if 
A<0. □ 

We first apply the Lemma to the case a — and a/3 > 0. It follows that the cto: cti and CT2 are uniquely determined 
if and only if — 1 < c and — 1 < d < 1. Therefore, the conditions —l<c,\d\<l are necessary for the well posedness 
of the initial-boundary value problem. This justifies the restriction on the matrix S = (Sj^^) to satisfy S'^S < 1. 
Next, we impose the Kreiss determinant condition [T^ ITsL IT^ which is stronger and requires that the constants ctq 
and (Ti can be bounded by the boundary data with a bound that is independent on w and C. It follows from Lemma 
^that this condition is satisfied if and only if |d| < 1. This immediately implies that the problem is well posed in 
if the initial data is trivial, since one can estimate the integration constants tro, cri by the boundary data (see Lemma 
8.4.3. in [3) while as = if -1 < c. 

On the other hand, suppose that d = and a — a which corresponds in setting to zero the transversal ingoing 
characteristic fields. The resulting initial-boundary value problem has been integrated in Ref. ITlj | by numerical means 
for different values of the parameters a and (3, subject to the restriction af3 > (in the notation of [Tll |. a = —71/2, 
P -272). In this case, d-C\~ d+{C^ - a) = -[(A + C)^ - 2a- l]/2 and since the map C ^ (^ + C)^ maps Re(C) > 
onto the outside of the unit disk minus the negative real axis, the system is ill posed if 2a + 1 > 1, that is, if a > 0. 
Indeed, the numerical results of Ref. 11] exhibit instabilities in those cases. 



IV. ENERGY ESTIMATES 



After the preliminary investigations in the previous section, we return to the case where e is a bounded 
domain with C°° boundary dfl, and derive some a priori estimates in this section. For simplicity, we only consider 
the source-free case here, where p = 0, Ji = 0; the generalization to the inhomogeneous case is discussed at the 
end of next section. We assume that we are given a smooth function u = {Ai,Ei, Wij) which satisfies the evolution 
equations 1)21314(1 . and the boundary conditions H14I15|I . The goal of this section is to show that this solution satisfies 
the estimates of Theorem |21 These estimates are then used in the next section in order to prove existence with 
semigroup methods. 

We start with estimates for the gauge invariant quantities Ei and Fij = Wij — Wji (i.e. the electric and magnetic 

fields), and the constraint variables C, Ck, C^i^\ and then estimate the gauge dependent quantities Ai and the 
symmetric part of Wij . It is the estimate for the latter quantities that requires a specific gauge condition, since the 
counterexample presented in sectionllllshows that the simple condition = does not allow to estimate the symmetric 
part of Wij in L/^ . Assume first that a = /3 = 0. In this case, the evolution system (|2I3I4|I is only weakly hyperbolic. 
However, in this case, one can replace Wij by Fij by taking the antisymmetric part of Eq. and discarding its 
symmetric part. Consider the boundary conditions H15|l : 

w^^ = 'w^^ + , w^'^ = 'E, ± n'^Fu , (59) 

where we recall that S^'' is a smooth matrix- valued function on dVt with the properties that S.^^'n^ — 0, S^'' Uj — 
and ^^j'^ — if 3i = and H^^ ^Sj^ < 5H'^^ for a 5 < 1 otherwise. The physical energy is defined by 

Taking a time derivative, using Eqs. H3I4|I and Gauss' theorem, we obtain 

= / {E^V'F,j+F'^V,E,)d^x= [ n'E^F,,dS 
Jn Jan 

= ll^^H^^{w'^^^w^+^~wt^w^r^)dS, 
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where dS denotes the surface element on 9r2. If the boundary conditions are homogeneous = 0) it follows 
immediately that the boundary i nteg ral is negative or zero, and we have an estimate for Ephys- If the boundary data 
is nonzero, we have the estimate [iJl 

H'^ (wf^w^^^ - w^r'^wY'^) < -eH'^ (wf^w'^Y^ + w^'^wY^^ + KW^g.g^ , (61) 

for some constants e = e{S) > 0, K = K{S) > from which we conclude 



for some constant a — a{S) > 0. 
When a(3 > 0, we obtain 



d 
Jt 



£phys{t) < EphysiO) + a I ||g(s)||i2(gn)C^S (62) 



Sphys = I {E^iV'F,,+aCj)+F'^V,E,)d^x 
Jn 

< a [ E^Cjd^x + a [ H'^g.^g^dS 
Jn Jan 



which contains the additional term J^^ Cj d'^x which cannot be estimated by Sphys alone. However, from the 
evolution system for the constraints, Eqs. (|7I8I9|I . and the constraint-preserving boundary condition 114II . 

y(+)=cy(-), y(±)^CT^n'=Cfc, (63) 

P 

where c is a smooth function on with |c| < 1, we can estimate the norm 

Scons (^niC('^)^^cf ' + T^C^ + C^Cfc) d^x. (64) 

Here, r > is a fixed constant with dimension of length and ni is a positive constant. We obtain 



/ h^'C^Y'^Cd^x-aT^ I n'^CCkdS 
2 Jn Jdn 



J, *-^cons 

at 



where we have used Schwarz' inequality and the inequality {h'-^ C^^^Y — 3C'^^'^^^ C'^^^ in the last step, and where we 
have set ni — 4a/ (3/3) and r = \J a.fi. On the other hand, using Schwarz' inequality again, we can estimate 



'J^E^C.d'x <^l^ f^E^E.+r'^C^C,^ d'x. 



Therefore, we obtain an estimate for the energy norm 

S — Sphys ~t" Scons ■> (^^) 

namely. 



d_ 

dt' 



£ <-£ + a [ H'^g^gj dS, (66) 
T Jdn 



which shows that we can estimate the gauge- invariant quantities Ej, Fij, C, Ck, C'^\ Notice that we can choose t 
arbitrary large, and thus we can make the exponential growth rate in the bound as small as we like. What is missing 
are estimates for the magnetic potential Aj and the symmetric part of Wij. Such estimates depend on the gauge 
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choice for (p. We have seen that the problem is not well posed in if we choose the temporal gauge = 0. An 
estimate can be obtained if we impose the following gauge choice instead: 

/^(jy^-V^Ek, n^V k(p ^ -n^ Ek on dO., (67) 

where A = V'^Vfe. The boundary condition on implies that on dVl, dt{n''Ak) = rJ^Ek + n^'Vk<t> = 0. This implies 
that = provided the initial data satisfies this condition. Next, we notice that we can estimate W since 

dtW = -(1 + 3/3/2)C+ A0 = -3/?C. Since we control C^^^ = W^j - V^Aj we can also estimate V*Aj. On the other 
hand, using Gauss' theorem twice, 

{V,A'){VjA^)d^x = / {VjA'){V,A^)d^x+ [ n' {A.VjA^ - A^VjA,) dS 
Ja Jdn 

3„ 



(yjA'){\'^A^)d'^x+ / K^jA'A^dS, (68) 
Jan 

where in the last term, Kij — {hik — ninkY^^rij denotes the extrinsic curvature of the boundary 9f2, and where we 
have used the fact that A^m = on dVl. Assuming that is positive semi-definite at each point of dVl^^, we obtain 
the inequality 



j (2V[M^1 • V[iv4j] + ViA' ■ d^x > j VM^ • ViA^- d^x 



(69) 



which allows us to estimate all spatial derivatives of Aj ; and in particular the symmetric part of Wij since we control 

Finally, we show how to estimate Af. Because of Eq. H67|l we have 

/ [E' + 2V'<t)){E,+2V,4>)d^x = i {E'E,+AE'V^<t) + AV''(t)V^4>]d^x 
Jn Jn 

= 4 / (j)n'{E,+ V,(/)) dS+ I {E'E, - Acj) {S/'E, + Acf) } d^x 
Jdn Jn 

E'E, d^x, (70) 



so we can estimate Ei + 2Vi0 and thus also dtAi — Ei + Vicj) since we have an estimate for Ei. The net result is the 
a priori estimate 



\Ht)\\lHn)+AHt)\\lHn)<^e'"^ 



\um\iHn)+r'\\vm\h(n)+ / llgWIlL^Oo)^^ 



(71) 



where a > 0, 6 > are two constants and, u = (r ^ Ai, Ej, Wij), u = (r ^C^J^"*, C, Cfc). This implies the estimate in 
Theorem 12 in the absence of source functions. 



V. EXISTENCE 



In this section, we prove well posedness for the initial-boundary value problem defined by the evolution equations 
(|1I2I3I4|) , the gauge condition p7|) and the constraint-preserving boundary conditions (|14I15II . The idea is to represent 
the problem as an abstract Cauchy problem 

^u{t) = Au{t), u{0) = uoe H, 

where A : D{A) C H ^ H is a linear operator on an appropriate Hilbert space H. This operator is basically given by 
the right-hand side of the evolution equations, and we will define its domain D{A) to be the space of smooth functions 
satisfying the boundary conditions with homogeneous boundary data. We then show that A has a unique extension 
A and that this extension generates a strongly continuous semigroup, which wc write formally as P{t) — exp(t^). 
Given initial data uq in the domain D{A) of the extension, the solution to the abstract Cauchy problem is given by 
u(t) = P{t)uo. The semigroup properties imply that ||M(i)|| < aexp(6t)||uo|| for some constants a > 0, 6 > which are 
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independent of mq; thus the problem is well posed. In the following, we assume that the coupling functions c and 
that appear in the boundary conditions are time- independent, lie in C°°{drt) and satisfy |c| < 1, W^S^^Sj'' < H^'^ 
at each point of dfl. 

The Hilbert space is motivated by the energy estimate in the previous section. Let 

Ha = {A e L'^{n,R^) : The derivatives VjAi exist in the weak 
sense and belong to i^(f2), and on dV,, n^Ai = 0.}, 
He = {'Ee L^(n,R^) : The derivative V'-Et exists in the weak 

sense and belongs to i^(ri).}, 
Hw = {W e L'^{n,M.^) : The derivatives VjW - V'Wji exist in the weak 

sense and belong to L^{il).}. 

Here, C — —V^Ei exists in the weak sense and belongs to i^(r2) means that C £ L^(il,R) has the property that 

n Jn 

for all test functions (j) £ C|J°(fi,R). Similar definitions apply for the weak derivatives VjAi and VjW — V^Wji. We 
introduce the following scalar products on Ha, He and H^: 

(A , B)a = / {r'^A^Bj + V'A^ ■ V,Bj) d^x, 
Jn 

[ {E^Fj + T^\/'E, ■ V^Fj) d^x, 
Jn 

{W,Y)w = I {W'^V,,+t''{V'W ^V,W'^){V,V ^V''V,k))d''x, 
Jn 

where A,B e Ha, E,F e He, and W, V e Hw- Notice that the requirement nM,; = on dVl makes sense because 
of the trace theorems. The following results are standard. 

Lemma 2 The spaces {Ha, {■ , ■)a), {He, (• , .)£;) and {Hw, (• , are Hilbert spaces. 

Lemma 3 Denote by C°°{Cl,W^) the class of functions 0, ~* R™ which are the restriction of a smooth function 

C°°(R3,R'") on n. Then, 

C°°{n,R) is dense m L^{n,R). 

{A e_C°°(f],R3) : n'A, = on dfl} is dense m Ha- 

C°°(fi,R3) is dense m He- 

C°°(rj,R9) is dense m Hw ■ 

We define the total Hilbert space H — i^(r2,R) x Ha x He x Hw with scalar product 

a 



(E, F)e - 



:,u)h = ^ (^T^p^ + T-^A'A^ + E'E^ + 2W^'^^W[ 



where u = {p, Ai, Ei,W,j), = Wij-V,Aj, C = p-V^E^, Ck = VkW-^^Wkj- The estimate ^ implies that 

the norm induced by (. , .)h is equivalent to the one induced by the scalar product constructed from (. , .)e^, {■ , ■)a, 
(. , .)e and (. , .)w Also, the Hilbert space H is topologically equivalent to the space H' defined in section^] Before 
we define the operator A on the Hilbert space H, we need the following result from elliptic theory (see, for example 

My- 

Lemma 4 Let F E C°°{ft), and g E C°°{dil). Then, the Neumann problem 

Au^F on n, (73) 
n'^VkU = g on dfl, (74) 

has a solution u E C°°{^1) if and only if 

'' Fd^x^ I gdS. (75) 
n Jon 

The solution is unique up to an additive constant. 
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This lemma allows to solve the boundary value problem (|67|l since according to Gauss' theorem 'V'^Ek (fix = 
Jg^ n^EkdS. We write (j) — —l\~^V^Ek + const for its solution. We define the linear operator A : D{A) C H H 
on the Hilbert space H by 

DiA) = {(p, A,E,W) e C°°{n,R^^) : On dn we haven' A, = 0, 



'E, + a,. 'n'^iWkj - W,k) - 0, W'^Ek -p + 6o- n'(V,W^ - 
where a = (a/) H(H - S)-i(H + S) and bo = (1 + c)/(l- c), and 



0}, 



(76) 



A 



(P 

A, 
E, 



\ 



Ej ~ \/jA-'^V'Ei 



(77) 



VWy - (1 + a)^SJ''Wj^ + aVjVF 

The existence proof relies on the following three propositions, which imply by the Lumer-Phillips theorem [T^ that 
A is closable and that its closure is the generator of a strongly continuous semigroup. 



Proposition 1 D{A) is dense in H . 



Proposition 2 The operator A : D{A) H defined in J 76| 77| ) is quasi- dissipative. That is, there is a constant b 
such that 

Re{u , Au)h < — {u , u)h (78) 
r 

for all ue D{A). 

Proposition 3 (A — A){D{A)) is dense in H for A > sufficiently large. 

Theorem 3 (Well posedness in the homogeneous case) The linear operator A : D{A) H is closable and its 
closure A is the generator of a strongly continuous semigroup P{t) in H. Given initial data uq G D{A), the map 
— > D{A), t P{t)uo is strongly differentiable and satisfies 



dt 



u{t) = Au{t), 



t > 0, 



(79) 



BiKE 



and so gives a solution to the constrained evolution system 
data, gi = 0. This solution obeys the estimate \iy\) . 

We will prove later that this solution also obeys the estimate (|2U|I . 



i|5l6)) with Jj — and homogeneous boundary 



A. Proof of Proposition ni 

The proof of Proposition ^ is based on Lemma 13 and 

Lemma h Let G e G°°{dn,m), F e C°°{dn,M?) with Fin* = 0, and let e > 0. There exists E £ G'^(n,R^) such 
that ||E||e < s and such that on dfl, 

H/E, = F, , V'^Ek = G, 

where Ui denotes the unit outward normal to CI, and — (5/ — UiU^ is the projection operator on the tangent space 
ofdCl. 

Proof: We construct E in the following way: First, extend ni to a neighborhood of dVl in by shooting geodesies 
through Ui at each point of dO, (such that n'^ViUk = 0). Denote by s the affine parameter which is such that V^s = ni 
and s = on dil. We consider a neighborhood Us of il which is spanned by s € {—S, 0) for some S > 0. Next, extend 
G and Fi inside this neighborhood by solving the ordinary differential equations 

V„G + kG' = 0, 

VnF,+KF,^ K,,F^, 

where the extrinsic curvature of the surfaces s — const. Notice that the second equation implies 

that V„(n*Fi) = —K{n^Fi), so n^Fi = on Us. Finally, let m > 2/6 and define the function e C°°(f2) by 
i^m = m~^ip(ms), where ip G C°°((— cx), 0], R) has the following properties: 
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(i) V(0) = 1, ip{s) = for all s < -2. 

(ii) ^'(0) = 1 

(iii) < -0'(s) < 1 for all s < 0. 

Then, we define Ei = Ui + e^rsV' ^^, where Ui — ipmniG, V — —'ipm^^'^^nkFi and £ijk denotes the natural volume 
element on il. By construction, we have 

so Ei satisfies the boundary conditions. Furthermore, 



lEll 



'Us 

< \Vol{U2/,n)\ 



^ / {G~V'F,fd^x+ f {F'F,+T^G^)d^ 



'Us JUs 

The right-hand side converges to zero as m — > oo, so the lemma follows 



□ 



Now let (p. A, E,W) e H, and let e > be arbitrarily small. According to Lemma|31 there exist p. A, E^^\W e 
C°°{n) with A„ = on dfl such that 



|A-A|U < 



E-EW|U<|, ||W-W||^<| 



5 ' " "-"5 

Using Lemma |31 we can find E(2) e C°°{n) such that HE'^'IU < e/5 and such that on dn, 



Therefore, (p,A,E(i) +E(2),W) e ^(yl) and 

Hp - p||i2 + ||A - A|U + ||E - e(i) - e(2)||£; + ||W - WIIh' < e. 
This proves Proposition ^ 

B. Proof of Proposition!^ 

Proposition El follows almost directly from the estimates in the previous section, so we only give the main steps 
here. Let u E D{A). Using Gauss' theorem and the estimate (|7U|) . we have, setting ip = —A^^W'^Ek, 

{u , Au)h = J^S^^A^iE, + 2V,<t>) + ^A^E, + aE^G, -\wG - '^h^^c'^^G^ d'x 

+ [ {2n'E^W[,j]-aT^n'GG^)dS 
Jan 



< 



ion 
1 

27 



n IT 



-A^Aj + E^Ej +r 



(3 ■' ZI3 4 3/? y 4 



where r = ^/af3, and Proposition |2] follows with 6 = (1 + r)/2. 
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C. Proof of Proposition 151 

Let V = (cr, B,F,V) e C°°{n,R'^^) with n'B, = on dQ, and let A > 0. We show that there exists u = 
(p, A, E, W) e D{A) such that (A — A)u — v; that is, such that 

Xp = a, (80) 
XAj = Ej - VjA-^V'^Ek + Bj , (81) 
XEj = V^VKy - Wj,) + a{VjW - V^Wji) + Fj , (82) 

XW^j - V.,E, + ^ %(V^£;fe - p) - V,V,A-iv'-i;fe + T/y . (83) 

A consequence of this is the eUiptic equation for E 

X^Ej = W,Ej + {af3 - l)VjV'E, + Sj , (84) 
where Sj — W^{Vij — Vji) + a{WjV — V*V,i) + XFj -— a(3X^^W ja, with the boundary conditions on dfl 

XH^ 'E, + a,'n\VkE, ~ V jEu) = -a^'n'^iVkj - V^k). 

r 



In appendix IXI it is proven in the case a^"* = a^H^^ where oq G C°°{dQ) is a strictly positive function, that 
for sufficiently large A > there exists a unique solution E e C°°(f2) to this problem. Setting p = X^^a and 
determining Aj and Wij from Eq. H81|) and Eq. H83|) . respectively, yields a solution m G E){A) to (A — y^)M = v. Since 
{(cr, B,F,V) e C°°(n,MiS) : n'Bi = on dn} is dense in H, proposition El follows. 



D. Intertwining operators and the constraint hypersurface 

Here we prove that the semigroup P{t) constructed in theorem |3| leaves the constraint manifold, defined by 

C = {u= [p, A„E„W,, ) e H : W^, = V,Aj , V'E, = p} 

invariant. In order to do so we introduce the Hilbert space G ^ {v ^ {c\Y\C,Ci) e L2(17,M13)} and the linear 
operator B : D{B) C G ^ G which is defined by 

D{B) = {(Cff \ C, Ck) e C°°(0,Mi3) : On dn we have C - bo^ n'C, = 0}, (85) 

P 

and 

-aV'^Ck I • (86) 

The operator B describes the evolution of the constraint variables, cf. Eqs. (|7I8I9|) . Similarly to before, one shows 
that D{B) is dense in G, that B is closable and that its closure generates a strongly continuous semigroup Q{t) in G. 
The Hilbert spaces H and G are related by the linear operator L : H ^ G which is defined by 

W,, - V,A, \ 

\ = \ P-^'E, . (87) 

Since L is continuous, the constraint hypersurface C = {u ^ H : L{u) = 0} is a closed subset of H . Furthermore, we 
have LD[A) C D{B), and we verify the intertwining relation 














H 


V J 





LAu = BLu, ueD{A). 



(88) 
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Since L is continuous and A and B are closable, it also follows that LD{A) C D{B) and that 

LAu = BLu, u£D{A). (89) 
This implies that for each u G D{A), the curve v : R+ D{B), 1 1-^ LP{t)u is difFerentiable and satisfies 

^v{t) ^ LAP{t)u = BLP{t)u = Bv{t), v{0) = Lu. 

Since B is the generator of the semigroup Q{t), it follows that LP{t)u — v{t) — Q{t)Lu, t > 0. Therefore, 

LP{t) = Qit)L, t > 0, (90) 
since D{A) is dense in H. Eq. (|90|l implies that P{t) leaves the constraint manifold C invariant. 

E. Proof of Theorem [H 

Finally, we discuss the case of nontrivial source functions and boundary data. We first reduce the problem to 
the case of homogeneous boundary conditions with nontrivial source functions and then use Duhamel's principle to 
construct the solution of the resulting problem. 

Let J(t) e C°°{n,m.^) cHE,t> O, be a continuous curve in iJ^, let uq = (p(°), A^, W^")) e C°°{n,R^'^)nH 
be smooth initial data, and let gi{t) e C°°{dil), i > 0, be a continuously differentiable curve in L^(9i7) with the 
property that gi{t)n'^ — on dO, for all t > Q. Assume that the initial and boundary data satisfy the zeroth order 
compatibility condition 

=Ef + a,^n'=(<) - wf^) = g„ (91) 



V'^^") - + n^vM"^ - "^'W^f) = 0, (92) 



on {t = 0} X dn. Choose ¥{t) e C~(0,R3) such that HjFj{t) = gjit) on dn and V''Fk{t) = on n (see Lemma 
O for all t > Q. We can choose F(t) such that it describes a continuously differentiable curve in _ff^(0,IR'^). Define 
v{t) = (0,0,Fi(t),0) e H, S{t) = (W,(t),0, J,(t),0) e H and J^(t) = Av{t) -v{t)+S{t), where Av{t) is defined by 
Eq. (|77() . One can verify that v{t), S{t), T{t), i > define continuous curves in the Hilbert space H . Next, consider 
the inhomogeneous Cauchy problem 

^w^Aw+T, (93) 
dt 

u;(0) = Wo = uq - w(0). (94) 

The compatibility conditions H91I92|I make sure that G D{A). We can solve this problem in an abstract way using 
Duhamel's principle: 

w{t) = P{t)wQ + I P{t- s)T{s)ds, (95) 
Jo 

where P(t) is the semigroup constructed in theorem and where the above integral is a TJ- valued Riemann integral. 
By construction, u{t) = v{t)+w{t) satisfies the inhomogeneous initial-boundary value problem (|1I2I3I4|I . ifTTjl . p4ll5|l 

with u(0) — uq: 

= Av - T + S + Aw + T ^ Au + S, (96) 

and u{t) satisfies the inhomogeneous boundary conditions. Furthermore, using the fact that Lv{t) = and LS{t) = 
for all t > 0, and using the intertwining relation 1)90(1 . we find Lu{t) = Q{t)Luo, which shows that the constraints are 
propagated. This concludes the proof of Theorem [21 
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VI. CONCLUSIONS 

In this article we considered the initial-boundary value problem for a formulation of Maxwell's equations which is 
structurally similar to the Einstein-Christoffel family of formulations of Einstein's field equations. In particular, our 
problem includes constraints that propagate with nontrivial speeds at timelike boundaries and so need to be controlled 
by specifying constraint-preserving boundary conditions. These conditions are supplemented with boundary conditions 
that control the incoming electromagnetic radiation. We have first considered the case of a fixed, a priori specified 
electromagnetic potential (f) and shown that the resulting initial-boundary value problem is not well posed in the 
expected sense, even though the Kreiss condition is verified. This is shown by an explicit counterexample representing 
an electrostatic solution in a "bad" gauge. This gauge is "bad" in the sense that it gives rise to modes growing linearly 
in Lut, where lu are the frequency components of the initial data. The counterexample also motivates how to deal 
with these modes by choosing a different gauge which eliminates them. In this new gauge, we are able to derive well 
posed initial-boundary value formulations; that is, we show that a solution exists in some appropriate Hilbert space, 
is unique, and depends continuously on the initial and boundary data. 

One key idea in our proof of well posedness is to start with an estimate for the physical energy instead of estimating 
the norm defined by the symmetrizer of the main evolution system, and with an estimate that controls the constraint 
variables. This gives rise to estimates for the gauge-invariant quantities. Once these quantities are under control, one 
estimates the gauge- dependent quantities (in our case the components of the vector potential and their first order 
spatial derivatives). These a priori estimates then motivate the Hilbert space in which the solutions are shown to 
exist using methods from semigroup theory. 

Unfortunately, the new gauge choice is of elliptic type, implying one has to solve a scalar elliptic equation at each 
time step of a numerical simulation. Although one can start an iterative method for solving it with the prior time 
value as seed, the procedure has some computational costs. Therefore, an interesting question is whether there exists 
other gauge choices leading to a well posed initial-boundary value problem which are easier to implement numerically. 
The relevant two places where the gauge choice is important is in estimating the norm of the vector potential 
Ai and of its derivatives ViAj. For the first case, since a estimate for Ei is given from the estimates for the 
gauge-invariant quantities, all what is needed is that the gradient of (j) be bounded in too. For the second case, 
the vanishing of the normal component. An, of Ai at the boundary is needed. In particular, this is true if An\dn = 
initially and n^{Vk<l) + Ek)\aQ, = for all times. So any gauge fixing satisfying these conditions should lead to a well 
posed system. 

How much of the above discussion can be carried over to General Relativity? As an example, consider the Einstein- 
Christoffel formulations with an a priori given shift vector that is tangential to the boundary. For this system, 
one has to specify six boundary conditions. There are three conditions that are needed in order to preserve the 
constraints Then, there are two conditions that should control incoming gravitational radiation. At least in the 
weak field approximation, where linearizations about Minkowski space are considered, it should be possible to specify 
such conditions in terms of the Weyl tensor, since - in the weak field approximation - this tensor is gauge invariant. 
Finally, there is a remaining boundary condition that has to be used in order to control a gauge freedom. We expect 
that for the a priori given shift case a similar counterexample as the one we discussed here can be found and that 
the resulting initial-boundary value formulation will not be well posed in the expected sense. On the other hand, 
we also expect that this problem can be avoided by requiring an appropriate elliptic equation for the shift and that 
well posedness can be derived for the resulting problem, at least for the case of linearizations about a stationary 
background. 

Finally, we briefly discuss a different toy model describing Maxwell's equations which resembles more closely the 
BSSN-type of formulations of General Relativity (see 21] for the original references and [H m |22b ,24j for an 
analysis of the mathematical structure of the equations). Instead of the variables Wij we introduce an extra variable 
r along with the constraint = F — V^A^ = 0, and consider the mixed first order-second order system 

dtA, = E,+V,^, (97) 
dtEj = W,Aj + aVjF -{1 + a)VjV'A, , (98) 
dtV = {l+l3)V''Ek+V''Vk(l), (99) 

where a, p are parameters subject to the condition a • /3 > 0. This model problem has been discussed by several 
authors ||^,2^,2^,27j in the past. The constraint variables C — —V^E^, propagate according to 



dtC = -aV'V.C^, 
dtC^ = -j3C. 



(100) 
(101) 
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In analogy to our previous model example, we consider the boundary conditions 

yW-cW-), y(±)^C'T^n'=VfeC^ (102) 
wi+^ = 'w^-'> , w^^'> = 'E, ± n'^iVkA, - V.^fe), (103) 
where |c| < 1, and S^'' is a smooth matrix-valued function on dVl with the properties that S^^'n^ — 0, S^^'rij — 



and H"^^ ^S^^ < H^". There is an analogous example to the one presented in section ITU which shows that it is not 
possible to prove well posedness in a space that controls the norm of the fields Ai, Ei, T and the first spatial 
derivatives of Ai when the gauge = is adopted. However, it is not difficult to see that one can adapt the estimates 
of section Hvl when the elliptic gauge condition H17|l is imposed instead, and so one can proceed as in this article to 
show well posedness. 
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APPENDIX A: ELLIPTIC PROBLEM 

Here we show existence of smooth solutions to the elliptic boundary value problem which arises in the proof of 
Proposition 121 We first simplify the problem and then apply standard methods to analyze the resulting boundary 
value problem. 

Let be a bounded open subset of M.^ with C°° boundary d^l. In the following, n = riidx* denotes the outward 
unit co-normal to d^}, and V denotes the covariant derivative with respect to the Eulerian metric ds^ — Sijdx^dx^ on 
M?" . We also introduce the projection operator Hi — Sf — riin^ on the tangent space of d^l, and the extrinsic curvature 
of the boundary surface, tHj = HI HjX/rUs. 

Let A and r be strictly positive constants, and let cq and d be a positive C°° function and a matrix-valued C°° 
function, respectively, defined on the boundary of fi, where d = (dj'^) satisfies n^d^^ = 0, d^^Uk — 0. We are 
interested in the following boundary value problem: 

X^Ej - WiEj - (r^ ~ l)\7jV'Ei = Sj on f7, (Al) 
{n^Vj + Xco)V'E, = G on dfl, (A2) 
n'{\/,Ej - ^jEi) + A dj'Ek = gk on dVl, (A3) 

where S € C°°(f],M3) and G € C°°(9ri,R), g e C°°{dn,m?) are given functions with the property that n'gj = 0. 
Notice that this has the form of the problem arising in Proposition O with r — Va^, 

1 1 1 -c 
rog r 1 + c 
d = = H(H + S)-i(H- S). 

In this appendix, we prove 

Theorem 4^ If cq > and d = c?oH where do is a strictly positive function on dfl, and if X > is sufficiently large, 
the boundary value problem I^A1\} . fAS^) . fASp possesses a unique solution E G C°°(f2,R'^). 
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We first simplify the problem and show that it is, basically, sufficient to consider the case where r = 1 and Cq is 
replaced by an arbitrary positive C°° function cp on 957. This freedom in the choice of cq will be important in view 
of constructing coercive quadratic forms. 

Let S e C°°(rj,R3) and G € C°°(ai7,K), g e C°°(arj,R3) be given. A solution to the problem (HTJ, (|X3|) 
can be obtained by performing the following steps: 

1. Decompose S = S^ + VJ where S"^ G C°°{n,M?), J e C°°(17,R) satisfy divS^ ee V*5f = and J\q^ = (this 
requires solving the Dirichlet problem AJ = divS, J|gQ = 0). 

2. Define G = (r^G + n^ViJ)\g^ and solve the problem 

(A^ - r^W,)(j) ^ J on 17, (A4) 
{n^Vj + A co)(t) = G on dfl, (A5) 

Using similar techniques as the one used below, it is not difficult to show that this problem has a unique solution 
<f> e C°°(17,M) if r > 0, Co > 0. 

3. Solve the boundary value problem 

X'^Fj - V^V.Fj = Sj on n, (A6) 

{n^Vj + A co)V'F, =0 on dn, (A7) 

n'(V,Fj -Vji^O + Ad/Ffe=gj-Ad/Vfc0 on .917, (A8) 

where cq is an arbitrary positive and smooth function on dil. Notice that a solution F S C°°(f2,R^) satisfies 
(A^ — V*Vi)divF = 0, {n^Vj + Aco)divF = Oj^^, so divF = according to the previous step. 

4. Then, it is easy to verify that E = F + Vi/) solves the boundary- value problem (|Aip . HA2p . (|A3I) . This 
solution is unique, because given two solutions their difference E'^ satisfies (A^ — r^V'Vi)divE'^ = 0, 
{n^Vj + Aco)divE'^ = 0|^^, so divE'^ = by the second step, and so E'^ satisfies the homogeneous version of 
the problem in the third step for which uniqueness will be shown. 

Therefore, it is sufficient to prove the statement of theorem 0] for the case where r ~ 1 and where cq is one particular 
positive and smooth function. Before we attack the problem, it turns out to be convenient to rewrite the boundary 
condition (jA2|) in a different form. Using Eq. (jAl|) with r = 1, we rewrite 

= [X'E^ - Sn) + W'^Vdn'VjEk - n^VkE^) + K'\V,Ek - VkE,) 
= [X^En - Sn) + W^V^igk -Xd^'Ej^), 

where En = E,n\ E]'^ = Hj'Ek, and where we have used the fact that is symmetric. Combining this with 
Ej = n^{n^V j + n)Ei + H'^^ViE^^ the boundary condition ljA2p can be rewritten in the form 

n\n'V, +n + Xc^')E, + {W'' - c^^d''')D,El' - c^^{Dd''')El' = Go , 

where Go — {X cq)~^ (G + Sn — D^gi) and where D denotes the covariant derivative with respect to the induced metric 
Hij on dfl. Finally, we combine this with the boundary condition ()A3|I and obtain the following boundary value 
problem: 

X'^Ej - V'VrEj ^ Sj on n, (A9) 
n^V.^Ej - DjEn + nj{W^ - c^^d^^)D,E*,f + m/'^^I-^ + vnjEn - c^'^nj{D,d'^)E[^ = n^Go + g^ on dVt, (AlO) 

where /ij^ = Xdjk + Hjk and rj — Acq"^ + k. 

From now on, we specialize to the case where d^j — doH^j where do is a strictly positive and smooth function on 
dn. Although we believe that this restriction is not necessary, it simplifies the proof below. Choose cq = do/ 2 in 
which case the boundary condition (|A10|) reduces to 

n'V.Ej - DjEn - njO'E*/ + jij^E{^ + 77 n^En - 2njd'^^{D'do)El^ = n^Go + g, on dVL, (All) 

where jMjk — Xd^Hji^ + and rj — 2XdQ^ + n. Notice that the third term in the boundary condition (jAllI) is minus 
the formal adjoint with respect to L^{dVt) of the second term in ljAll|) . This and the positivity of jiij and 77 will allow 
us to find the appropriate elliptic estimates and to prove 
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Theorem 5 Let do € C°°(9f2) be strictly positive at each point of dfl. Denote by ki < K2 the nonzero eigenvalues of 
Kij at each point of dfl, and set N = ^ H^i Di log do ■ Dj log dp . Then, if X is sufficiently large such that 

A > max{0, ^ [N - + «2)rfo) (A12) 
Z do 

the boundary value problem 1^A9\) . l^A 1 1]) possesses a unique solution E e C°°(f2,M"^). 

The proof uses standard methods [sOl [IJ : We first recast the boundary value problem in weak form and show 
the existence and uniqueness of weak solutions. After that, we derive regularity results which show that this solution 
is, in fact, smooth. 



1. Weak formulation 



In the following, let fc > 0, denote the Sobolev space of functions M™ whose partial derivatives 

of order < k are in L'^{Q). We abbreviate W'' — iJ*''(ri, R-^). Introduce the linear operators L : , 
b : ^ 7jfe+i/2 i^QQ^ defined by 

i(E), - iX'-V^\7,)E,, 

6(E)j = n'V.Ej - DjEr, - n^D^El^ + ti/E^f +r]njE„ - 2njd^^{D'do)El^ , 

We consider the boundary value problem i(E)j — Sj, fe(E)j — gj, where S € W'' and g e H''+^^'^{dfl,R^). 
Let E,F g C°°{n,R^). Using Gauss' theorem we find 

/ F^L{F,)jd^x = [ {X^F^Ej + \/'F^ ■ \/,Ej) d^x - [ F^n'\7,EjdS 
Jn Jn Jan 

= Qa(E,F)- / F^b{E),dS 
Jan 

where dS denotes the surface element on dQ and where we have defined the bilinear form 

Qa(E,F) = [ {X^E^Fj + \/'E^ ■S/iFj)d^x 
Jn 

+ [ [fi'^ E,Fj+rjEnFn-2FniD' log do)El^ -FnD' El' -F^'D' En] dS. 
Jan 

Therefore, 

(F, L(E))L2^n) + (F, 6(E))i2(9n) = Qa(E, F). (A13) 

Because of the trace theorems we can extend Qx to x W^: If E,F € , they are also in H^/'^{d^), and (in 
the sense of distributions) DiEn £ H~^^^{dfl), which is the dual space of so the integral Jg^F^'D'EndS 

makes sense. Furthermore, there is a constant C such that 

|gA(E,F)| < C||E|lvvi • \\F\\w^ for aU E,F e W\ (A14) 

Thus Qx is a bounded bilinear form on the Hilbert space W^. Equation (|Il3|l . which holds for aWE e W^, F e W^, 
implies 

Lemma 6 Let E g W'^, and suppose that S G , g g H^^'^{dfl,W^) are given. Then, F is a (strong) solution of 
n^) . I^A11\) if and only if 

(F, S)mn) + (F, g)man) = Qa(E, F). (A15) 

for all F £W\ 

Proof: The "only if" part is clear. In order to show the "if" part assume that l|Il5|) holds for all F g W ^. In 
particular, it holds for all F g Cg°{n,M.^). Since Cg°{n,M.^) is dense in L'^{n,M.^) it follows, in view of Eq. (|IT3|) 
that L(E) = S. Next, choose F g C°°{n,R^) with Fj ^ at the boundary, which implies b{E)j ^ g^. □ 

This motivates the following 
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Definition 1 Suppose A > 0, and let S e , g G R^) 6e ^wen. We call E e a weak solution of 

i/ equation l^Al^) holds for all ¥ eW^ . 

The next lemma implies the existence of weak solutions: 

Lemma 7 Suppose the assumptions of Theorem [3 are met. Then, Q\ is coercive. That is, there exists a constant 
0-0 = 0.0 (A) > such that 

Qa(E,E) >ao||E||2^,l (A16) 

for aWEe . 

Proof: Let E e C°°(f7,R3). We have 

gA(E,E) = [ {X^E^Ej+S/'E^ ■\/iEj)d^x 
Jn 

+ I [fi'^E.E, +riEl - 2E^{D' logdo)^f - E^D'E]^ - E'^D'E,,] dS. 
Jan 

The last two terms in the boundary integral form a total divergence, and so their contribution to the integral vanish. 
Therefore, setting ao = min{A^, 1}, = ^y Wi Di log do • Dj log do and using Schwarz' inequality, 

Qa(E,E) > ao||E||^i + / {fi'^E.E, + E^ ~ N [doH'^E.Ej + d^^El])dS 

Jon 

= ao||E||^i + / ([(A - N)doH'^ + k'^ E,E, + [(2A - iV)d^i + k] E^) dS. 
Jan 

> ao||E||^i. 

Since C°°(f2,M'^) is dense in and since Q\ is bounded, the lemma follows. □ 

Corollary 1 (Existence and Uniqueness of weak solutions) Suppose the assumptions of Theorem O are met. 
Then, the problem fA9\} . I^A 1 1\) possesses a unique weak solution. 

Proof: Define the linear functional J € W"^* by 

J(F) - (F, + (F, g)L2(aa) , ¥ 

According to the Lax-Milgram lemma there is a unique E such that Qa(E, .) ~ J. q 

For later use, we introduce the operator Cx : ~f W^* which is defined by 

(£aE)(F) = Qa(E,F) 

for E, F e W\ Since |Qa(E,F)| < C||E||vi/i • \\F\\w^ it follows that WCxEWw^- < C\\F,\\w^, so Cx is bounded. 
Furthermore, it is injective because of (£aE)(E) — Qa(E,E) > ao||E||^i and onto because of the implications of the 
Lax-Milgram lemma. The weak form of (jA9|) . |IA11|I is simply £aE — J. 
The following estimate will be important: 

Lemma 8 Cx ■ — > W^* defines a linear bounded and bijective operator. For each E G it satisfies 

Wniw^ < ao'llCxniw^' ■ (A17) 
Proof: According to Lemma we have, for E G W^, 

aolMw^ < Qa(E,E) = (/:aE)(E) < ||/:aE||v^i. • \\E\\w^ . 

□ 
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2. Regularity 

Theorem 6 (Regularity) Suppose the assumptions of Theorem\^are satisfied. Let k >0 and suppose that S G W'', 
g ^ H'^'^^^^ (dn,M.^) . Then, the unique weak solution E of the boundary value problem fA9j) . 1 1]) lies in W^^'^ . 
Furthermore, there is a constant Cfc > such that 

||E||vi/fc+2 < Ck (||S||vi/fc + ||g||/i-fc+i/2(aQ) + ||E||vyfc+i) . (A18) 

Remark: For fc = in particular, the theorem shows that any weak sohition hes in and so is a strong solution 
according to Lemma |S1 

Remark; Since the solution is unique, one can show (see, for example. Lemma 8.38 in Ref. |^) that the last term 
in the right-hand side of the estimate l|A18|l can be dropped. The resulting estimate then implies well posedness in 
the sense that the solution depends continuously on the data. 

We break the proof of the theorem into several steps. First, we choose a finite covering of and show that the 
problem can be localized. Applying difference quotients to the estimates (|A17|I . (|A18|I we then show interior regularity 
first. Next, we flatten the boundary by choosing an appropriate coordinate patch and map the problem to a half plane 
problem. Using difference quotients again, we first show differentiability in the directions tangential to the boundary, 
and then use the solution properties to show differentiability in the normal directions as well. 

We start with the case fc = 0. 



Localization 



Suppose E G is a weak solution, and let x be the restriction on of a function in Cg°(M'^). Then x • F € 
for each ¥ & and J(xF) = (/:aE)(xF) = (£a(xE))(F) +i?(E,F), where, using the Leibnitz rule and Gauss' 
theorem, one finds 

i?(E,F) = Qa(E,xF)-Qa(xE,F) 
Jn 

+ f [D'x{FnE'/ + Er,F'J)~WnX-E'F,]dS. 
Jaa 



Therefore, >Ca(xE) = J where 



with 



J(F)-(F,S)i2(a) + (F,g)i2(ao) 



S, = x^, - 2V,x • y'E, ~ iW^x)Ej , (Af 9) 

9j = X9j + [V„x • E, - D,x ■ En - n.D'x ■ E'J'] , (A20) 

Since E S W^, we have S € W'^ and g e i?^/^((9r2, R'^), and S, g have the same support as x- Therefore, it is sufficient 
to assume that E and the data F, g are supported in some coordinate patch. 
Finally, we notice that taking F g C,5"(r2,R'^) in equation (jAf 5p implies that 

(A^ - V'Vi)Ej = Sj (A2f ) 

in the sense of distributions. 



4. Interior regularity 

Assume first that E and S are supported in Q (and so the boundary data vanishes). Then, we can extend them to 
elements in iJi(R3,R3), and iJ"(R3,R3), respectively By rewriting (A^ - V*V,,)i^j = Sj in the form (f - W,)E^ = 
Sj + (1 - X'^)Ej, it follows immediately from the next lemma that E e H'^{R^,R^). 
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Lemma 9 For each s gM., the operator 1 — V'Vi induces an isometry of H''{R^) into H'^ ^(M'^). 

Proof: This follows immediately using Fom'ier transformation. |— I 

In view of the generalization to operators with variable coefRcient and in preparation of the boundary regularity 
result, we give a different proof which is based on the use of difference quotients. So let E S 7f^(M^,R'^) be a weak 
solution with compact support. The idea is to apply the estimate (IA17II to the difference quotient -Df E, defined by 

D'lEix) = i (E(a; + /le,) - E(a;)) , 

where h > Q and denote the unit standard vectors in . The partial derivative DiE, exists in the weak sense if 
and only if D^E is bounded in h [sj- According to Eq. (|A17p . we have 

\\D^E\\H.<a,'\\C^D'iE\\H., . (A22) 

Now for any F G iJi(M'\R3), 

{CxD'iE){E) = Qa(D?E, F) = Qa(E, Dv'^E) + i?(E, F), 

where 

R{E,E) ^ Qy,{D'lE,E) - Qx{E,Dr'^E). 

Using the "Leibnitz rule" and "integration by parts" for differences quotients it is not difficult to see that there 
is a constant C > (which depends only on the metric coefficients and their derivatives) such that |i?(E,F)| < 
C||E||//i • ||F||//i. Furthermore, since E is a weak solution, 

Qx{E,D-''E) - J{DT>'E) 

< ||S||l2(r3) • ||L),j^''F||i2(R3) 

< ||S|U2(R3) • ||F|lffi(a3) . (A23) 

Taking everything together, we obtain 

||I?fE|l^i(R3) < 1 (|lS||i2(R3) + C||E||^i(K3)) . (A24) 

Taking the limit /i ^ we obtain DiE E i/^(M'^,R'^). Repeating the argument for each i = 1,2,3, it follows that 
E € aj^j that 

||E||i/2(R3) < Cl (||S||i2(R3) + ||E||^1(R3)) , (A25) 

for some constant Ci. 



5. Boundary regularity 

Let Xo € drt and suppose that E and the data S, g vanish for all |x — xo| > i5 > 0. If 6 is sufficiently small, we can 
introduce new coordinates y such that the boundary is described by y = = (fiiV^iy^)- That is, we can 

consider the problem on the half plane M!^. In order to show regularity, we proceed as above, except that now we are 
only allowed to take the difference quotients D\E, A = 2, 3 in the directions tangent to the boundary. Also, instead 
of Eq. HA23p . we have 

Qx{E,D^''F) = J{D^''F) 

< I|S||l2(r.3 ) • ||i:)^''F||^2(R.3 ) + ||g||^l/2(R2) • ||-D^''F||h-1/2(R2) 

< (||S||l2(r3_) + ||g||Hl/2(R2)^ l|F||ffl(R3_) , 

where we have used the trace theorem in the last step. Therefore, it follows that DaE e i?^(Mi]_, R"^), A = 2,3. In 
order to show that DiE g i/^(R']_, R'^) we remember that L{E) = S in the weak (distributional) sense which allows 
us to express (£'i)^E in terms of all other second (and lower order) derivatives of E for which we know that they are 
in L'^{R\). Therefore, we have E e i72(R3^,R3), with an estimate 

||E||^2(r3_) < Cl (||S||^2(r3_) + ||g||^l/2(R2) + | ] E | ] ^1 (r3_ ) ^ . (A26) 

Taking a finite covering of il with a corresponding partition of unity, and summing up the interior and boundary 
regularity results, it follows that E e W'^, and the estimate (jA18|) holds for A: = 0. This proves theorem for A: = 0. 
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6. Higher order regularity 

Suppose E G W^+^, fc > 1, is a strong solution of L(E)j = S,, b(Ey = , where S g g e H''+^^^{dn,R^). 
Suppose also that we have the estimate 

||E||^y.+i < Ck-i (llL(E)ll^.-i + ||5,(E)|| llElli^.) (A27) 

for all E e W'^^^ . We show that this implies that E e W''^'^. The statement of theorem |S1 then follows by induction. 

The idea of the proof is the same as for k — 0: We first localize the problem: Let x be the restriction on of 
a function in C^{R^). Then x • E g W'''^'^ satisfies the problem L{xE)j = Sj, b{xEy = g' where S S and 
g g H'^+^/'^idU.^R?) are given by Eqs. ljA19|) and (|A20p . respectively. So we can assume that E and the data are 
supported in a coordinate patch, and so it is sufficient to consider the problem on R'^ or on the half plane M^. 

We discuss only the case Ri^ here; interior higher regularity follows in the same way. We apply the estimate (|A27|I 
to D^E, where A = 2,3 denote tangential directions. Introducing the commutators 

i?''(E) = ip^E) - L>^L(E), 
r''(E) = h{D'\^)- D'Xhi-E.), 

which satisfy the estimates (which follow from the "Leibnitz rule" for difference quotients) 

\\R^{^)\\H^-Hm.\) < C'I|E||//)c+i(r3 ) , 

lk''(E)||_H-fc-l/2(R2) < C||E||^fc + l/2(R2) , 

where the constant C depends on bounds for the metric coefficients hij and their fc + 2 order derivatives, we obtain 
ll-DAE||_H-fc+i(iB;3) < Ck (^||Z)^L(E)||^fc-i(jj,3 ) + ||-D^&(E)j||^fc-i/2(R2) + ||E||^fc+i(R3^)^ 

< Ck (^||S||//fc(R3_) + ||g||^fc + l/2(R2-) + ||E||//fc + l(R3 )^ . 

Therefore, DaE g iJ'=+i(R^, M^)^ A = 2,3. Because E satisfies the equation i(E) = S this also impHes DiE g 
i?''+i(R3^,IR3)^ so E g iJ'=+2(R3 ^R3) ^j^jj corresponding estimate. 
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